hep-th/0305018 

The Fall of Stringy de Sitter 

Andrew R. FreyQ Matthew Lippert^ and Brook William^] 
Department of Physics 
University of California 
Santa Barbara, CA 93106, USA 
(Dated: February 1, 2008) 

Abstract 

Kachru, Kallosh, Linde, & Trivedi recently constructed a four-dimensional de Sitter compact- 
ification of IIB string theory, which they showed to be metastable in agreement with general 
arguments about de Sitter spacetimes in quantum gravity. In this paper, we describe how discrete 
flux choices lead to a closely-spaced set of vacua and explore various decay channels. We find that 
in many situations NS5-brane meditated decays which exchange NSNS 3-form flux for D3-branes 
are comparatively very fast. 
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I. INSTABILITIES OF DS FLUX COMPACTIFICATIONS 



De Sitter spacetime (dS) holds a special place in the study of quantum gravity. Construct- 
ing and exploring the maximally symmetric spacetime with positive cosmological constant 
A has been the source of much recent interest despite (or perhaps because of) its stubborn 
opacity. While much progress has been made in the understanding of the other maximally 
symmetric solutions, Minkowski and anti-de Sitter (AdS) spaces, dS has until recently eluded 
string theoretic description because of some of its unique properties. The observer-dependent 
horizon of dS, like a black hole horizon, yields a thermal state with finite entropy. Not only 
are the S-matrix observables of string theory precluded in A > spaces 111], but, due to the 

□ n 

inevitable Poincare recurrences |2j, all observables are ill-defined |3|]. These issues would 
merely be of abstract theoretical importance were it not for recent observational evidence 
jj] indicating that not only was A > during in the early universe during inflation, but it 
seems to be so today. 

It has become increasingly clear that dS cannot be a stable state in any theory of quantum 
gravity. The symmetries of dS are incommensurate with the discrete spectrum implied by 
finite entropy |5[ . Rather than a stable vacuum, dS is instead a metastable resonance whose 
lifetime, on general entropic grounds, must be less than the recurrence time P, 0, lzj| . 

One would ask, then, what does string theory say about dS and its decay modes? String 
models of dS have been difficult to find partly because, as non-supersymmetric vacua, they 

moduli stabilized. Notably, some dS compactifi- 
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are isolated points in moduli space with al 

cations of string theory were described in |8. 0| and. in a woll-controllod manner lor critical 
strings, by Kachru, Kallosh, Linde, & Trivedi (KKLT) |6(. Generically, any string theoretic 
dS compactification can decay and decompactify because the 10D Poincare invariant 

string vacuum is supersymmetric and so has vanishing energy density. However, this is far 
from the only decay mode. For example, in any compactification in which RR fluxes con- 
tribute to the potential, D-brane instantons change the fluxes and the cosmological constant. 
This has been an object of study in many papers, including [ill . 

We consider a slight twist on the brane instanton decays. In p], the cosmological con- 
stant gets a positive contribution from D3-branes, rather than directly from the fluxes. This 
effect has been seen in the AdS/CFT correspondence, where instantonic NS5-branes provide 
a decay mode for the D3-branes [r|. These results apply to the similar dS compactifica- 
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tions of KKLT and are particularly of interest because they can end in a state of positive 
cosmological constant. Therefore, one might wonder whether this type of decay could occur 
quickly enough to affect the cosmological constant within the age of the universe. In this 
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to dS compactifi- 



paper, we generalize the results of Kachru, Pearson, & Verlinde (KPV) 
cations and compare the decay rate through the 5-brane channel to two other decays, one to 
decompactification and the other by D3-brane tunneling in the compactification manifold. 
We give explicit examples in which the 5-brane decays are much faster than the others. 

In the next section, we review the dS vacuum construction that we will study. In section 
TTT1 we flesh out the discrete landscape of vacua that are available through tuning and among 
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and 



which our instantons will interpolate. We then review the AdS/CFT instantons of 
make the corrections necessary to compactify their backgrounds in section IIVI We apply 
our calculation to find decay times for specific sets of initial parameters in section and 
compare them to those of KKLT in section IVII In addition, we comment on two other 
possible decay channels. We will generally keep factors of the gravitational coupling k 4 
and the string length a' explicit in formulae, but any numbers we cite should be taken in 
Planck/string units. 

II. BUILDING DS VACUA 

Constructing a solution of string or M-theory with a four-dimensional dS vacuum has 
been a longstanding challenge. Such a solution must be non-supersymmetric and requires 
aspects of the theory beyond the low-energy SUGRA limit. 

Recently, however, KKLT j3] presented a specific construction in critical string theory 
with no unfixed moduli. The model was based on the warped flux compactifications studied 
by Giddings, Kachru, & Polchinski (GKP) [17] 1 . Non-perturbative corrections fix the overall 
Kahler modulus of this tree-level no-scale model, resulting in a stable, supersymmetric AdS 
vacuum. KKLT then added D3-branes to yield a metastable dS vacuum and showed, by 
considering decays to decompactification, the lifetime to be less than the Poincare recurrence 
time. 

1 The GKP type of compactification was studied earlier in simpler cases and in M-theory by lsl llfli EH 
2.1] . The supersymmetry conditions and equations of motion were considered in [^H EH Y2A \. Explicit 
constructions on tori and K3 are in 25, 26, 271. 
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The GKP compactification of IIB string theory on a threefold M with 7-branes and 03- 
planes can be efficiently described as an F-theory compactification on a CY fourfold (CY) 
X. X is elliptically fibered over M such that the fiber's complex structure r = cq + ie~^ is 
the IIB axion-dilaton (we take for simplicity r = i/g s ). We will consider the orientifold limit 
of F-theory in which M is an orientifolded CY threefold. Three-form fluxes and D3-branes 
are added subject to the global tadpole constraint, or the global conservation of RR 5-form 
F 5 flux: 

= - N m + — I H 3 AF 3 - ^— ^ . (1) 
2k? /x 3 J M 3 3 24 V > 

The Euler number of the CY fourfold x{X) gives the effective negative D3-brane charge in 
IIB of 03-planes and D 7-branes wrapped on 4-cycles of M. For typical choices of X, xPO 
can be up to O(10 5 



28| . This is must be balanced by the charge from 4D space- filling 
D3-branes, D3-branes, and the wrapped NSNS and RR 3- form fluxes H 3 and F 3 , which also 
source F 5 . 

To construct their model, KKLT began by choosing X and a set of wrapped fluxes, while 
setting A^D3 = = 0. The CY threefold M has 63 ^> 1 three-cycles, and a particular choice 
fluxes if 3, F 3 G H 3 (M, Z) represents a point in a 263 dimensional lattice. The fluxes combine 
into a single complex 3-form G 3 = F 3 — tH 3 . For simplicity, KKLT chose h 1,l (X) = 2, so 
that M has a single Kahler modulus p. In addition to the moduli r and p, M has h 2,1 (M) 
complex structure moduli z a . 

In the presence of fluxes, the classical 4D effective Af = 1 superpotential is 

w = ^ [ G 3 A n , (2) 

K 4 J M 

where Q is the holomorphic (3,0) form on M. Wo then is given by the (0,3) part of the G 3 
flux which, because the fluxes are quantized, can only be tuned discretely. The tree-level 
Kahler potential (ignoring warping of the spacetime metric) 

K = -3 logH(p - p)) - log(-z(r - f )) - log (-\ f O A Cl) (3) 

along with Wq gives the no-scale potential 

V = e K WDiWDjW (4) 

where i,j sum over all moduli but p, fCy = didjK, is the Kahler metric, and Di = di + diK is 
the Kahler covariant derivative. Except for the volume modulus p, this potential generically 
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fixes all other moduli such that G3 is imaginary self-dual. 2 The remaining condition for su- 
persymmetry, D p W = is satisfied only when Wq = 0, which implies that in supersymmetric 
vacua G3 is a (2,1) form. 

The geometry of M is, of course, very complicated but is accurately described near 



conifold points by the Klebanov-Strassler (KS) solution 
deform the conifold; at the tip y — 0, the metric is 



. Wrapped fluxes warp and 



ds 2 = h- 1 /\ v dx^dx v + bg s Ma' {e 2u dy 2 + dft 2 + e 2u y 2 dQ 2 2 ) (5) 

where b ~ 1 is a numerical constant and e u is the compactification length scale (here we use 
10D string frame). Notice that the S 3 at the tip has a fixed proper size depending only on 
the fluxes. Also, the S 2 is nontrivially fibered over the S 3 . Away from the tip, the throat 
has approximately a warped conifold metric 

ds 2 « hr^^dx^dx" + h 1/2 e 2u (dr 2 + r 2 ds 2 Thl ) (6) 

where ds^ M is the metric on the base T 1 ' 1 . In this region, the warp factor is approximately 

h = l + (L 4 /r 4 ) log(r/r s ) (7) 

with the length scale L 4 = |e" 4tt g s M«'. Here, the radial coordinates r and y are compli- 
cated functions of each other, and the tip is at y — 0, r = f. For the undeformed conifold, 
the singular tip is located at r = r s = fe _1//4 . Splitting the conifold into the tip and throat 
in this manner is described in Q m d refa enee s therein. 

The radial modulus Imp = e 4u /g s = o is defined so that, at large radius, the total 



unwarped volume of the compactification is J M d 6 x\/h l / 2 g ~ e 6u a' 3 . The fluxes through 
any 3-cycle of M are quantized, and for a given conifold throat 

M =7^J f ^3, K = -^- [ H 3 (8) 



B 



Air 2 a' J a ' Aix 2 a' 

where the A cycle is the S 3 which stays finite at the tip and the B cycle is the six-dimensional 
dual of A. GKP found that the warp factor at the tip of the conifold is related to the 



2 It is natural to wonder if corrections to the Kahler potential due to warping could fix the radial modulus. 
While the precise form of JC is difficult to compute for p, because the 10D solution exists at tree level 
for all compactification scales, the final potential must be no-scale |3 We will look at warping in the 
complex structure Kahler potential below. 
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deformation parameter z of the tip, which is determined by the flux superpotential (0), by 
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h(y = 0) « f^ 3 , * = exp 



2tcK 



It is this particular form of the warp factor that gives the A cycle a fixed proper size at the 
tip. Note that this is not the r — > r limit of equation ((7j) because the conifold is deformed. 



To generate a nontrivial potential for p, as suggested in |17(, KKLT considered non- 
perturbative corrections to the superpotential (J2J). Both wrapped Euclidean D3-branes and 
gluino condensation on the worldvolume of non-Abelian D7-branes generate additional terms 
of the form 

5W = Ae iap (10) 

where the constants A ~ 0(1) and a ~ O(10 _1 ). For simplicity, KKLT took p to be purely 
imaginary, p = ia, and A, a, W to be real. The potential now becomes 

V = ^^{Ae^(l + ^) + W }, (11) 

and for suitable Wq < there is a supersymmetric vacuum with Vo < 0, implying the non- 
compact directions are AdS. For \Wq\ -C 1, the AdS minimum lies at a cr ^> 1 where the 
SUGRA can be trusted and a' corrections are small. 

The final step in the KKLT construction is to add enough D3-branes so that Vq > 
and the vacuum is dS. The global F 5 charge must still be conserved via eqn ((H), and the 
addition of p D3-branes gives = —p. By adjusting the fluxes, a corresponding increase 



in J M H 3 A F 3 balances this reduction. The D3-branes break supersymmetry and add some 

m 

SV=^; D = 2p, 3 h-\r) (12) 
a 6 



IM 

extra energy 



where p 3 is the brane charge. To minimize their energy, the D3-branes migrate to a coni- 
fold tip, so the energy density per D3-brane depends, through eqn @, on the fluxes. For 
sufficiently fine-tuned parameters, this additional term in the potential lifts the AdS global 
minimum to a dS local minimum. 

Unlike the AdS vacuum, the dS minimum is only metastable. KKLT investigated one 
possible decay mode, tunneling to large a. The potential becomes arbitrarily close to zero 
at large radius, so it is possible to tunnel to a runaway, decompactifying solution. 

Coleman and De Luccia (CDL) [22| described such an instanton including gravitational 
back-reaction. In terms of a canonical scalar field (p = (a/3/2 log a) / /t 4 , the Euclidean action 
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S E [<p]= I '^Xy/g(~R + ^(dip) 2 +V{ip)) 

J 4 

= - J d 4 x^gV{<p) (13) 

using Einstein's equations to get the second line. The instanton (Pcdl is an 0(4)-symmetric 
interpolation between the dS vacuum at <p cr and the supersymmetric vacuum at ip = oo. 
When Wick rotated back, this gives the usual expanding bubble of true vacuum inside the 
false dS vacuum. The action of the static dS vacuum is simply computed to give 



So = = So (14) 

where Sq is the entropy of the dS vacuum. The tunneling probability per unit volume 
is given by the difference between the action of the instanton solution and the static dS 
vacuum: 

P%% ~ e^™** . (15) 

From eqn (fT3|) . S[ip] < for V(<p) > 0, and the resulting lifetime is exponentially less than 
the Poincare recurrence time t r ~ e °: 

^~C*-«<tr (16) 

which is in line with the general arguments of 

In addition to the CDL instanton, KKLT considered decompactification decay via the 



stochastic Hawking-Moss (_HM) instanton [34J. Considering decays of general dS string 
compactifications, [5] and |7( also discussed thermal fluctuations using the HM reasoning. 
Whereas the CDL instanton tunnels through the potential barrier, the HM instanton relies 
on thermal fluctuations to carry ip to the top of the potential, where it can then roll down 
the other side to the true vacuum. While the original HM process is homogeneous, KKLT 
argued it should be interpreted as a horizon-sized fluctuation. If the potential has a broad, 
flat maximum at tp\, the state there is approximately dS with energy V(ipi) > Vq and entropy 
S±. The probability per unit volume for a thermal fluctuation is given by the difference in 
entropies between the fluctuation and equilibrium: 

P^ay - e S *- S ° . (17) 
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The decay time t^ ay = (P^ay)' 1 * s again l ess than the recurrence time t r and is also less 
than the CDL decay time (jlfij) when the potential barrier is short and wide and thus the 
thin-wall approximation is invalid. 



III. FINDING DS PARAMETERS 

As described in section |Hj obtaining the vacua constructed in [6j requires fine tuning 
subject to several constraints. First, one must adjust the bulk fluxes so that \Wq\ <C 1. 
Moreover, a dS minimum requires fine tuning of the fluxes, K and M, in the KS throat. A 
given value of Wq tightly constrains one's choice for Dp (cf. eqn IT2*|) . For example, KKLT 
presented a model with Wq = — 10 -4 and an AdS minimum of Vq = —2.00 x 10~ 15 ; by 
adding one D3-brane with D = 3 x 10 -9 , they achieved a dS minimum of Vq = 1.77 x 10~ 17 . 
This is a very special choice of fluxes indeed. For Dp < 3 x 10~ 9 the minimum is at 
Vo < and is AdS, and for Dp > 7.5 x 10~ 9 a local minimum no longer exists. There are 
additional constraints as well. In 16] it is shown that there exists a classical instability if 
p/M > 0.08. Furthermore, results from section 3 of ^| rely on approximations valid when 
K/{g s M) > 1/2. 

With such fine tuning and taking into account that the tuning parameters K and M are 
discrete, one might question if it is possible to build such a model at all. Such tuning would 
require the existence of a "discretuum" 3 . We have done numerical searches in order to map 
out the discrete landscape of dS vacua. Figure [T] shows the existence of the discretuum. 
Here we have plotted the possible values Vq that have a dS minimum and can be achieved 
with integer fluxes for the parameters used in KKLT, Wq = — 10 -4 , a = .1, A = 1, g s = 
0.1, K4 = a' = l. 4 It is clear that for a desired value of Vq there exists a configuration of 
fluxes with Vq = Vq + e, where e is very small; i.e. a discretuum does exist! For each of the 
models studied K/g s M > 1/2. Here we have allowed M to range from 75 to 1000. The lower 
bound avoids the classical instability (for p < 6) As one goes to higher and higher values 
of M, one must also increase the amount of induced D3-brane charge on the D7-branes in 



3 The authors of 01 coined this term to refer to situations in which a discrete spectrum is sufficiently dense 
to allow for an (almost) arbitrarily fine tuning. Our discretuum is not as finely spaced as those in 

4 In addition to tuning Vq by varying the fluxes M and K, one could, in principle, vary Wo by adjusting the 
bulk fluxes. While this would certainly increase the discretuum density, we leave Wq constant as explicit 
calculation of Wq in terms of bulk fluxes is prohibitively complicated. 




FIG. 1: The possible dS vacua with Vq for given M illustrate the density of states consistent with 
a discretuum. 

order to satisfy This might require adding more D7-branes and, thus, more degrees of 
freedom, which, though massive, could cause problems when considering loop corrections. 

The smallest possible value of Vq, for the parameters used in KKLT, is O(10~ 20 ), a far 
cry from the desired O(10 -120 ). In order to obtain a more realistic vacuum energy, one must 
attempt to construct a background with |Wq| ~ O(10~ 55 ). While such a fine tuning seems 
improbable, with 6 3 sufficiently large, it is at least possible 5 , if not particularly natural [13]. 

We have so far considered only a single KS throat, as in KKLT. However, a general CY 
has many of them. By considering backgrounds with multiple KS throats the discretuum 
density is increased dramatically. One finds that (|12|) becomes, 

tfv = y-Aft. Di = 2 ^h-\f l ). (18) 

i 

where % labels the different throats. Clearly, by adjusting the fluxes in each individual 
throat, one may tune SV with greater accuracy. For a single KS throat we found O(10 3 ) 
configurations with a dS minimum. Analogously, for 2 KS throats (75 < M\ < M2, 75 < 
M 2 < 300) we find O(10 5 ) dS minima. It is easy to find configurations with O(10) KS 
throats 6 , leading to an amazingly dense set of vacua. The inclusion of a second throat also 

5 One can estimate the smallest \Wo\ to have log(|Wo|) ~ —263. We thank S. Kachru for discussion on this 
point. 

6 For example, in |35l | a family of quintics are constructed with 16 conifold singularities. 
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not worry about these subtleties (other than the well-known effects on the instanton 
32I ]) since our main focus is on instanton decays ending in dS. The configuration with 



lowers our minimum value of Vq by an order of magnitude. Though this is nice, it does little 
good in helping build a model with a realistic cosmo logical constant. We suspect that even 
with the addition of 10 or more throats the lofty goal of Vo ~ 10~ 120 would still be far out 
of reach. 

The following sections describe various decays, analogous to those studied in [r| , in which 
one unit of ^3 flux is exchanged for M D3-branes. For geometries with single KS throats, 
after one decay the final state has a negative cosmological constant and a big crunch in its 
future. It has been argued that these decays should not be allowed in a quantum theory of 
gravity and also that instantons mediating these decays may not be possible to construct [36] . 
We wil 
action 

multiple KS throats is more interesting. As with the single throat, these may decay directly 
into states with negative cosmological constant. However, there can now be decays from 
one dS vacuum to another with smaller A (modulo some classical evolution we will discuss 
later). This process is of particular interest, since it allows for a rather generic set of fluxes 
on several KS throats to undergo a series of decays to dS vacua with smaller and smaller 
cosmological constant; this situation is similar to that envisioned by 0] and expanded upon 

by HQQ- 

IV. DECAYS A LA KPV 

A. Review of NS5-brane Instantons 

The KS geometry found at conifold points of GKP compactifications was first studied in 
the usual decoupling limit of string-gauge theory dualities j^. The relevant gauge theory 
dual is a duality cascade with an energy dependent effective number of D3-branes; in the 
IR, most of the D3-branes have been transformed into 3-form fluxes. The BPS domain wall 
that transforms the D3-branes to fluxes was described by KPV; it is a polarized NS5-brane 
that carries D3-brane charges and bends over the A cycle at the deformed conifold tip [16]. 
As the NS5-brane moves over the A cycle, the D3-branes are absorbed into the background 
RR flux, and the background NSNS flux jumps by a unit due to the NS5-brane charge. 

KPV also described nonsupersymmetric gauge theories with p D3-branes at the tip of the 
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FIG. 2: Top: In the KPV process, p D3-branes polarize into an NS5-brane wrapping an S 2 on the 
A cycle. The NS5-brane then slides to the opposite pole, becoming M — p D3-branes. Bottom: 
In the thin- wall limit, the NS5-brane is instead wraps the A cycle at a particular Euclidean radius. 

conifold, as in KKLT. Due to the 3-form flux background, the D3-branes suffer a classical 
instability to brane polarization (first discovered in J38I]) as an NS5-brane wrapping an S 2 in 
the A cycle. However, for p > M/12, the NS5-brane itself is unstable to collapse around the 
A cycle, reducing the NSNS flux and turning the D3-branes into supersymmetric D3-branes. 
For smaller p, the decay of the NS5-brane proceeds by tunneling; in Euclidean spacetime, 
the NS5-brane is slightly polarized with D3 charge at infinity and bends around the A cycle 
to leave D3-branes at the origin Q]. This process is illustrated in the top line of figure El 

For p small enough, KPV showed numerically that the thin-wall approximation is very 
reasonable. In that limit, the instanton appears to be an NS5-brane wrapping the full S 3 
of the A cycle at a fixed radius, as shown in the bottom line of figure 121 The wrapped F 3 
flux induces M units of charge in the NS5 worldvolume gauge theory which is canceled by 
the charge carried by the ends of M D3-branes. The p D3-branes end on the outside of the 
NS5-brane, and M — p D3-branes end on the inside. The bubble tension in the effective 
theory is just the NS5-brane tension times the volume of the A cycle. These instantons are 
clearly related to the BPS domain walls KPV found. 

In the rest of this paper, we will focus on the thin-wall limit to estimate the instanton 
bubble tension. As has been argued strenuously [36|, the thin- wall limit certainly does not 
describe the full picture of the decays, but the other contributions to the Euclidean path 
integral (such as Hawking-Moss instantons at the other extreme) should only enhance the 
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decay rate. Therefore, we take the point of view that the thin-wall limit estimates an upper 
limit for the decay time. As a consequence of the thin-wall limit, we may, as in KKLT, 
ignore the polarization of the D3-branes in the initial metastable state. Before we turn to 
the modifications necessary for including KPV instantons in KKLT compactifications, let 
us also note that our instantons are cousins of the supersymmetry-changing domain wall 
bubbles found in [3], just as the AdS/CFT instantons of KPV are related to the BPS 
domain walls. 



B. Corrections for Compactifications 

There are several modifications that we have to make to the KPV instanton decay formula 
due to the fact that we have a compact GKP geometry rather than a noncompact conifold. 
The first and most obvious correction is that gravity is no longer decoupled, so we should 



include the effects of gravitation on the decay time. These effects are well known 
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12. 



321 ] : in appendix El we work out the specific formula we need. The decay time, including 
gravity (but ignoring the large number of massive fields in the compactification), is tdecay ~ 
exp[— ASe], where ASe is the difference of the Euclidean actions for the instanton and the 
initial background state as given in eqn (jB10|) . It depends only on the bubble tension, the 
initial vacuum energy density, and the change in energy density. Given two dS states from 
section ITTTJ we just need to calculate the bubble tension and plug into (jB10|) . 

There are also modifications to the tension of the bubble. The easiest to calculate is an 
effect of working in the 4D Einstein frame. Let us emphasize that we need to work in the 
4D Einstein frame to use the superpotential formalism of section |Hj and this is also the 
frame in which the potential has been calculated. The Einstein frame is also the frame used 
in calculating the instanton decay time. It is easiest to get this by going to the NS5-brane 
action 

Se = -f / d^x^/det g^S I d 3 x^/g^ , (19) 
9s J J 

where g^ is the 4D pullback of the 10D metric, 5 is the delta function at the radius of 
the bubble (with the determinant of the metric included), and #53 is the determinant of 
the metric on the A cycle. The 10D string frame and 4D Einstein frame are related by 
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h l l 2 g^ u = g s 2 e 6u g fiU , so the NS5 action becomes 

/ 7 2/3 \ 3 / 2 1,3/2 _ 

8. = r^^[— ) (2^)(^Ma-)^= ifa3a , 3/2gyV/4 . (20) 

(Henceforth r is the instanton bubble tension.) In the first equality for the tension r 5 , we 
have separated the contribution from the conversion to Einstein frame, the warp factor, and 
the volume of the A cycle. We have ignored the contribution to the action from the NSNS 
6-form potential, which KPV showed is negligible in the thin-wall limit. Heuristically this 
is because the 6-form potential only has two legs on the A cycle and the 5-brane fills the 
entire cycle, as shown in figure El However, the RR field strength F 3 gives a worldvolume 
anomaly that requires M D3-branes to attach to the 5-brane. Here, there are p D3-branes 
on the outside and M — p D3-branes on the inside. 

The other correction we should make is due to the action for the moduli. Since the moduli 
are fixed by the flux superpotential (J2J), after the NS5-brane bubble changes the flux, the 
VEVs of the moduli will change. Therefore, we need to take into account the rolling of the 
moduli to the new vacuum. We will focus on the deformation modulus z of the conifold 
for the following reasons. First, it clearly changes significantly when K changes (see Q). 
Also, for a noncompact conifold, K does not affect the dilaton or other moduli, so we would 
expect that they would be only minimally affected by a change of K in the compact case 
(the other moduli are typically fixed by fluxes on other cycles). Also, KKLT have shown 
that the VEV of a does not change much due to the presence of D3-branes. Therefore, 
since we expect g s and a to keep roughly the same values before and after the decay, we 
expect that they will not roll much, and we will treat them as constants. There is actually a 
significant tree-level potential for g s and a when z is not at its VEV, and we will consider its 
effects in the next subsection. Nevertheless, we expect our estimate of the contribution from 
z not to be affected significantly by other moduli. To be conservative, one could multiply 
the contribution from z by a fudge factor, but we note that we are only making an estimate 
to begin with, so we are not quite that careful. 

To estimate the tension due to the rolling of z, we will assume that just inside the NS5- 
brane z is in its original vacuum value outside of the NS5-brane and rolls quickly to the new 
VEV inside. This is probably not the exact classical solution, but we will use it and the 
thin-wall approximation as an upper limit. At tree level (where we are working), we can 
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write the action as 



Se(z) = \ f d*x y /gi[]C tS d ll zd ,l z + K\e K K. zS D z WD 2 W] 

271-2 r 

= — d£r 3 [K, z ~ z (K eg dsz - nle K/2 - luJ D s W) (K gg d^z - n 2 A e K / 2+iu> D Z W) 

+K 2 4 e' c/2+iuJ d i :zD z W + Kle^-^d^zDsW] , (21) 



where u is some phase (physically, we have to take it so that the Euclidean action comes 
out positive because it started positive definite). As above, r is the radius of curvature of 
the bubble, while £ is the radial coordinate corresponding to proper distance. This is clearly 
minimized when only the last two terms contribute. Taking the average Kahler potential in 
the exponential, we get (up to numerical factors of order unity) 

r z ^27r 2 e {lc)/2 (\AW\ + \A)C\\(W)\) . (22) 

(This comes from the definition of the covariant derivative and the chain rule.) This deriva- 
tion is very similar to that of BPS domain walls and is also used in j^j]. Actually, it is 
easy to generalize this estimate to include other moduli, but we will only consider z in the 
superpotential and Kahler potential. We should note that AW and A/C are calculated from 
the inside of the NS5-brane (where z is not in a vacuum state) to the new vacuum on the 
interior of the instanton and not from the original vacuum to the new vacuum. Since we are 
just making an estimate, (• • • ) will be an average value over the region of variation of z. 

The change in the superpotential is given entirely by the superpotential of the conifold 
just inside the NS5-brane minus Wq. This is because in the vacuum states, the K and M 
fluxes are (2,1) forms and so do not contribute to the superpotential (see, for example, j^ll). 
Using the notation and conventions of 12, Q' we S e ^ 

AW = -W{z) = ^\^1JJ± ( MQ{z) - i—z 

~ W°* : f» fa -i)j"»'fV (23) 



«| \27Ti g s J g s Kl 

where K is the NSNS flux on the inside of the bubble and z is evaluated outside the bubble. 
This follows from the definitions 

[ fl = a ' 3/2 z , / O = a' 3/2 Q(z) , Q = — zlogz + holomorphic (24) 

J A JB 27T2 
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and the relation from eqn (JSJ) that z(outside) ~ exp(—2ir / g s M)z(inside) . To overestimate 
(W), we will take 

« |W | + |AW| . (25) 

The Kahler potential is significantly more complicated, and, because we are concerned 
with a modulus that lives at the bottom of a throat, we need to take the warp factor into 
account. Including warping and bunching the Kahler potential for all other complex moduli 
together into K c (that is, integrals over other cycles), eqn (JHJ) becomes 



K,(complex) = — log 



K>4 \Ja J B J A JB 



(26) 



To compute j B Qh, we use the trick that the cycles have a monodromy B — > B + A around 
z = (in the same way it was used to find the leading term in Q) and eqn Q (which is 
valid at points both inside and outside the bubble) to expand out 

JC(complex) &K C - e ^ a '^ 9s ^ |z| 2 / 3 log \z\ 2 . (27) 

(This is, to our knowledge, the first calculation of part of a Kahler potential with warping 
included.) Actually, there will be other terms in the B cycle integral, but it is reasonable 
to believe that, as in the unwarped case, this is the leading term that depends on z. Then, 
using ® and assuming the complex structure gives small contributions to /C c , we get roughly 



16(T 3 ' 



AJC 



a'\g s Mf 2/3 
2txk\ 1 1 



Air 
~9sM 



D 4Tr/3g 3 M 



+ log|z| 2 (e 4 ^ M -l) 



(28) 



The total bubble tension is therefore 



T 5 + T z 



b^ 2 z 2n 2 gl /2 



+ 



167r 3 a' 3 / 2 ^ 5/ V 9 / 4 4crV 2 \ 9s ^ 



(2tt)V 5 / 2 



z + (W + 



(2vr) 2 « 



2^/5/2 



9s K 4 



X 



a'\g s Mf /3 



47T 



3 47r/3g s M 



g s M 



+ log \z\ 2 (e 



in/3g a M 



1) 



(29) 



With the bubble tension in hand we are now in a position to calculate decay rates. However, 
before moving on we would like to take a closer look at subtle issues ignored in the above 
calculation. The anxious reader, fretting over the fate of his or her universe, may skip ahead 
to section |3 and leave the following subsection for a more careful reading. 
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C. Other Considerations 



1. D3-brane Migration 

The KPV instanton bubble not only reduces the NSNS flux K and annihilates D3-branes, 
but it also leaves behind D3-branes. If there are D3-branes in other throats, the D3-branes 
will feel an attraction and roll through the bulk 7 and into the throat with the D3-branes. 
Eventually, they will annihilate with the D3-branes via tachyon condensation. If this mi- 
gration is part of the instanton, then, in many cases, all the D3-branes will be annihilated, 
leaving a Big Crunch spacetime with negative energy density. If there are more D3-branes 
to start, the final state could still be dS. 

However, we argue that we should not consider the migration of the D3-branes to be part 
of the instanton, but rather as a classical process that occurs after the bubble nucleates. Our 



logic is something like the discussion of the bounce instanton of quantum mechanics in 42] ; 
the instanton should only tunnel through the barrier to some energy slightly lower than the 
initial state, and classical evolution should take over. Typically, in the thin-wall limit, we 
just assume that the inside of the instanton is just the final state. In our case, though, we 
expect that the D3 migration would not be well approximated at all by a thin-wall instanton 
because they are very far from the D3-branes, so the potential is very flat. (Contrast this to 
the case for the z modulus, where the gradient of the potential is Planck scale.) This logic 



is consistent with the discussion of Hawking-Moss and related instantons in 

We expect the migration time to be similar to the bubble thickness for the motion of 
the D3-branes in the Euclidean description of the instanton. The migration times are larger 
than the bubble radius for the rest of the instanton, so we will treat the D3-brane migration 
as a classical process. In fact, the migration times are larger than the initial dS radius itself 
for the models we consider, which is the maximum bubble radius. 

In appendix El we estimate the classical migration time for a single D3-brane migrating 
from one tip to another. For the particular model we examine, AtM ~ O(10 15 ) (in string 
units). As discussed below, decay times for the instantons we are considering are much 
larger, O(exp[10 9 ]). Thus, inspite of the fact that total migration will vary a great deal 



7 We assume forces due to objects in the bulk, such as D7-branes with gluino condensation [41|, can be 
ignored. We thank S. Kachru and L. McAllister for discussion on this point. 
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from model to model, the total decay time At TOT = At rfecaj/ + At M ~ At decay is relatively 
unaffected. 

We should note that the classical D3-brane migration followed by D3/D3 annihilation 
could leave a state with negative cosmological constant. In cosmology, if the spatial slices 
have nonnegative curvature, the FRW constraint equation means that the universe cannot 
actually transition to a negative cosmological constant. Instead, there is a Big Crunch 
singularity j3, Q • Though it is preferable to end in a dS state after the full decay, 
this is not necessary as long as the initial instanton has a lifetime much longer that the age 
of the universe. Note that the instanton, however, ends in a state of positive cosmological 
constant, so we avoid the concerns raised by 



2. Rolling Radius 

Now we should go back and examine the classical potential for a that arises because z 
is away from its VEV. The behavior of the radial modulus in flux-generated potentials has 
been studied in an attempt to find inflationary behavior in we are in a different regime 
here because we do not take z to be slowly rolling. One point to address is that we cannot 
actually calculate the Kahler potential with warping for z excited because it is not clear if 
eqn @ would still hold as z changes. However, we will assume that it is valid since the 
starting and ending points of our evolution are vacuum states for some values of the flux K. 
The key point is that for instantons that go from dS to dS, the boundary conditions on a 
mean it should not roll much, so the following discussion does not apply. What we are doing 
here is comparing instantons with different boundary conditions, one with a unchanged in 
the final state and one with a — > oo in the final state. 

We make the comparison as follows. The classical potential for a and z naturally pushes 
a to large radius as long as z is not in its vacuum state (note that this potential is extremely 
large compared to the KKLT potential (jlljl . so we can ignore the KKLT potential here). 
We will make a very rough estimate of the change in a while z rolls to its vacuum. If we 
believe that a changes enough to get over the barrier of the KKLT potential before z reaches 
its vacuum and the classical potential vanishes, then we expect dS to Minkowski decays - 
mediated by NS5-branes! - will dominate over dS to dS decays. This is because the classical 
evolution should have a lower action. Otherwise, the dS to dS decays will dominate, at 
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least in the NS5-brane channel. We will not say anything else about these dS to Minkowski 
decays since they are less computationally tractable and are somewhat redundant with other 
decays to large radius. 

Now we can roughly estimate the potential for a and z. As in GKP [17], we work assuming 
small z, which implies that d z W > (d z JC)W, so we will consider only the derivative of the 
superpotential. As before, the D P W terms cancel with — 3|1U| 2 . As a final approximation, 
we take only the leading terms of the Kahler metric for z small. Thus, we approximate the 
potential as 

. .. n/2 lrl 4 /3 M K 2 

(30) 



V = gte- 12u (2n) 5 
We have used 



a' 2 


z\ 


4/3 


4(g s M) 2 


log 


\z 
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M . K 2 

— \ogz + — 
2vr g s 



(g s M) 2 a' 3 4/3 

as the Kahler metric for z. This is singular at z — 0, but our evolution never takes z — > 0. 

To get a very rough estimate of the change in radial modulus u (remember that a = 
e Au /g s ) while z changes, we approximate that the proper distance in the u direction of 
moduli space is proportional to the proper distance moved in the z direction of moduli space. 
The proportionality constant is given by the directional derivative (in the moduli space 
orthonormal frame) of the potential. Using the Kahler metric (|31jl to get the orthonormal 
frame, we find that 

Am w — — 1 =^- (32) 



Au w ( ^ sM) Y% | 2 / 3 | log|z| 2 | (e 2 ^ M - 1) (33) 

187T Kj^ 



up to factors of order unity. We have used y/lC^Az for the proper distance in the z direction. 
The factor of \/l2 in (|3*2*|l comes from the normalization of u. 

Using the potential graphed in KKLT as a guide, we expect that Au only needs to be 
> 0.1 for the Minkowski decay to predominate, which is achieved by z > 10~ 3 . As it turns 
out, we will mainly be interested in cases with smaller z, so we will not consider the 5-brane 
mediated dS to Minkowski decays any further. 
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3. Thermal Enhancements 



Due to the fact that dS has a temperature, we might expect that the 5-branes that make 
up our instantons should have some nonzero entropy. Since the exponential of the entropy 
gives a density of states, the decay time should be reduced by a factor exp[— S(NS5)]. This 
argument was first given in Q|. There it was argued that the brane instantons probably 
are out of thermal equilibrium with any matter or radiation in the cosmology, so they 
should have a temperature corresponding to the dS temperature. However, whether the 
temperature should be the initial dS temperature, final dS temperature, or the geometric 



mean was undetermined. It is now clear 



471 ] that the brane would be in equilibrium with 



the initial dS because it corresponds to accelerating observers in the two dS spacetimes. 8 

We will, however, neglect this effect. The bubble temperature is just the inverse radius, 
T = W Q. The^e, the tem pe rature is no. M g n enough to excite the "Kaluza- 
Klein" modes of the bubble much, and the entropy would access only the zero-mode quantum 
mechanics. We expect that the enhancement factor would be relatively weak, therefore. 

V. CALCULATION OF DECAY TIMES 

Throughout this paper, we have mainly discussed the KPV instantons as CDL thin-wall 
instantons. However, they contain an NS5-brane, which makes them also of the membrane 



class of instantons studied by 



11 



12j . In appendix El we demonstrate the equivalence of 
these two formalisms by showing that they give the same decay rate given initial and final 
cosmological constants and instanton tension. 

Using the results of section HYI and appendix El we are able to calculate decay rates. 
For illustrative purposes let's first consider a model with a single KS throat. In particular, 
for 3 D3-branes sitting at the tip of a throat with K = 12, M — 87, one finds that the 
probability per unit volume for NS5-brane mediated decay is P ~ exp(— 10 19 ). Decays to 
decompactification are much faster, P ~ exp(— 10 17 ). We expect this to generally be the 
case for single throat models. Moreover, as discussed in 36], since all single throat decays 
will have A < in the final state, the instantons mediating these decays might not exist. 
It is for this reason we have chosen to focus on models with 2 KS throats, which, after the 



We thank the authors of |4?J for sharing their results with us prior to publication. 
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Model 


Pl,P2 


Ki 


Mi 


K 2 


M 2 


z\ x 10 17 


z 2 x 10 5 


AA x 10 31 A_ 


. x 10 17 


1 


1,1 


9 


15 


3 


19 


4.2 


4.9 


3.9 


69 


2 


1,1 


9 


15 


4 


26 


4.2 


6.3 


3.9 


2.7 


3 


1,1 


9 


15 


9 


69 


4.2 


28 


3.9 


4.5 


4 


1,5 


9 


15 


8 


51 


4.2 


5.2 


3.9 


4.5 


5 


1,5 


9 


15 


13 


91 


4.2 


13 


3.9 


7.6 



TABLE I: Models and Cosmological Constants 



Model 


(t/t c ) 




-9 


KKLT: (t/t c ) 


ln(Sf ) >< I " 18 


T2T: 


ir/r c ) 


Ht T d Zy) X ID" 18 


1 


0.163 


0.66 




1.8 


0.32 




24970 


0.35 


2 


0.164 


86 




7.7 


8.9 




24257 


8.9 


3 


0.164 


40 




5.9 


5.2 




16512 


5.2 


4 


0.163 


3.7 




2.9 


1.1 




34054 


1.1 


5 


0.164 


18 




4.6 


3.1 




33517 


3.1 



TABLE II: Tensions and Decay Times 

initial decay, have A > 0. 

What follows is a discussion of the decay rates for several different two throat models. 
Table |l] shows the the fluxes and number of D3-branes, Pi, in each throat. In each model the 
initial KPV instanton occurs in throat 1. This decay is driven by the notably small value of 
z±, which makes the tension very small. Note that we have specifically chosen models where 
this is the case. The change in and resulting value of the cosmological constant (AA and A_ 
respectively), due to KPV decay, are also given in table E The small value of z\ corresponds 
to small |AA|, which would increase the decay time, but this effect is compensated by the 
small bubble tension. How the decay rate depends on these values is given explicitly by eqn 
(IB10|) . We list the tensions and decay times 9 for the KPV instantons in tabled along with 
tensions and decay times for two other decay modes discussed in section IVII below. Note 
that the lifetimes for these models are ~ exp(10 9 ), where as the age of the universe (times 

9 Note that these are the decay times for a unit volume, i.e. t decay = f -1 
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the horizon volume) is ~ exp(10 3 ), so even the most anxious reader can now relax and enjoy 
the rest of the paper. 

For each of the models discussed above, although the initial instanton decay yields a 
spacetime with positive cosmo logical constant, the ensuing D3-brane migration results in a 
negative cosmological constant, a situation which, as discussed in ultimately 
leads to a Big Crunch singularity. Note, as previously mentioned, this is a classical process 
and thus avoids arguments given against instanton decays to negative A j^. The total 
migration time, as shown in appendix El is negligible compared to the decay time 10 and will 
thus be ignored. We should also note that, although it seems difficult to find two throat 
models with a positive cosmological constant after D3/D3 annihilation, it should be possible 
to construct multiple (> 2) throat models that end in dS. 

VI. COMPARISON TO OTHER DECAY MODES 

The KPV instanton is just one of several avenues by which these dS vacua can decay. 
One particular mode, thoroughly studied in (f| Q| and reviewed at the end of section HU is 
tunneling to decompactification (in the CDL formalism). In these decays, or for any decay 
in which A_ = 0, ASe takes a particularly simple form, 

ASe = ~ n 2,2 ■ ( 34 ) 

(1 + t-2/t 2 ) 

For comparison purposes, we have calculated the CDL tensions and decay times (t^^ T ) 
for five models discussed above. These are also listed in table |H] Note that in each model 
the tensions are super-critical, r/r c > 1. This will in fact always be true for decays to 
decompactification since, 

t 1 r°° 

- = ...... WQ / dtP^ZVW) > 1 , (35) 

r c v/4V(0+)/3 

for any V(<fi) whose barrier width (in string/Planck units for our normalization) is greater 



that a/2/3. Noting that So < 0, it is clear from (J33|) that the lifetime, t^ ecay ~ exp(— ASe), 
increases with r/r c . Though the story is more complicated when comparing to decays with 
A_ 7^ 0, this will generally still be the case, and it is this fact which drives tfJ^ T to be much 



10 Note, however, that it is long compared to the string scale, O(10 15 
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greater that t*fj^. Take careful note that table [XT] lists the logs of the decay times. For these 
models t demy T / ^ decay ~ ex P (10 8 )! These KPV instantons are, in technical terms, much much 
much faster. It is possible to find super-critical KPV instantons in which t^ e f a y T < t^^f . 
However, these require larger z in the decaying throat, leading to larger initial cosmo logical 
constant and slower decay times. 

Another particularly simple decay mode occurs in models with multiple KS throats. The 
potential energy of a D3-brane is proportional to ft, -1 (r), the inverse warp factor given by eqn 
(JHJ), which is locally minimized at the tip of each throat. However, the energy is lower still 
at the tip of other throats with smaller z. D3-branes can therefore tunnel from one throat 
to another. On the other hand, h ~ 1 in the bulk, presenting a substantial potential barrier 
;hrough which to tunnel. These instantons are similar to the glueball decays considered in 

mw\. 

As in previous examples, we consider models with two KS throats. The D3-brane portion 
of the total potential is initially (cf. eqn (fTSjO 

SV = 2 -^h-\h) Pi + 2j ^h-\f 2 ) P2 ■ (36) 
a 6 a 6 

After the tunneling occurs, the form of 8V is unchanged except for p\ — > p\ + 1 and p 2 — > 
p 2 — 1. These decays have little effect on a, and thus o will be treated as a constant 
throughout this calculation. 

To find the decay rate, we compute the instanton tension from the Euclidean brane action 
in the thin- wall limit using j^: 



T = ^ttv^P^/V jT dr 2^^[h-\r)-h-\h)] ■ (37) 

The prefactor is from the conversion between r and a canonically normalized scalar in the 4D 
Einstein frame. Note that here we are using rescaled coordinates so e 2u does not appear in 
the metric ©• We then plug r and A± = n 2 A (V + 5V±) into eqn ()B10|) to obtain the "throat- 
to-throat" decay time {t^ 2 J ay ). Once again, the these decay times and tensions are listed 
in table |H] Note, however, that these instantons tunnel to negative cosmological constant; 

.„ r „„&.. „„,«„ „„«,„.,„,,.. 

of [32j. As with the KKLT decays, the tensions are uniformly supercritical, and t^ ay is 
remarkably similar to tf e f a y T . Indeed, we expect these decays to be super-critical because 
the conifold throats are long in string units, giving a wide potential barrier. Moreover, 

22 



taking the limit n\r ^> A + in eqn ()B10|) . it is easy to see that, 

24tt 2 

(AS E ) T2T » — - , (38) 
A + 

and thus from (JUj) and one can see that £^; T ~ t™ a y 

The reader should remember that we are working only in the thin-wall limit and that 
Hawking-Moss instantons can also give significant contributions to the decay rate. However, 
in the models we have described, the dS to dS decays are subcritical, so the Hawking-Moss 
contributions seem unlikely to change our qualitative results; KKLT found that Hawking- 
Moss instantons begin to dominate over thin-wall instantons only when r ~ yV{}Pi) > t c . 

The KPV instanton deals only with changes to fluxes and branes in eqn One might 
speculate about processes which could involve changes to x(X), or the induced D3 charge 
on wrapped (p, q) 7-branes in the IIB. From the F-theory viewpoint this would obviously 
involve topology change. While one could consider non-trivial D7-brane worldvolume gauge 
fields undergoing a small instanton transition and emitting D3-branes into the bulk, we 
know of no analog for non-trivial curvature on a four-cycle. A possible x changing instanton 
would involve D7-branes unwrapping a particular 4-cycle and wrapping a different one; 
however, these 4-cycles would be homologous unless the D7-brane can tear, so the induced 
D3-brane charge would remain the same. However, it may be interesting to explore whether 
X-changing mechanisms are possible. 



VII. THE END OF THE WORLD (AND THIS PAPER) AS WE KNOW IT 

The decays considered in this paper in a very real sense would represent the end of 
the universe for anyone unfortunate enough to experience one. 11 Note, however, unlike 
the decays in CDL, even when D3/D3 annihilation following a KPV decay results in a 
Big Crunch, lifeforms might be capable of knowing joy for 10 -28 seconds while the D3- 
branes migrate across the compact manifold. We can all take comfort in the fact that 
even the fastest decays we consider have decay times incredibly greater than the age of our 
universe. Assuming that our calculations hold even approximately for a compactification 
with a realistic cosmo logical constant, we will have to worry about the death of the sun long 
before the death of the universe. 

11 Apologies to REM for the section title. 
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Of interest, however, is the fact that we constructed decay modes other than the straight- 
M deca y to defalcation diseased in Q fl fl. In fact, we found it ea Sy to 
construct NS5-brane mediated decays that occur much more rapidly than the decompact- 
ification decays. We reiterate that the NS5-brane decays can have a subcritical tension. 
It is also noteworthy that the final state of many decays is not 10D Minkowski spacetime 
is instead dS or a space with negative cosmological constant which ends in a Big Crunch. 
In fact, depending on the region of parameter space, we found that decays mediated by 
NS5-branes can end in dS, 10D Minkowski, or with negative cosmological constant, without 
considering other decay channels! The lesson is that, even in the KKLT models, there are 
many different metastable vacua and many different possible decay modes. 
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APPENDIX A: MIGRATION OF D3-BRANES 



The KPV decay leaves D3-branes at the tip of the KS throat where it occurs. In this 
appendix we analyze their subsequent classical motion in configurations with multiple KS 
throats. The D3-branes, produced by a decay in one throat (Throat 1), are attracted 
by D3-branes in another throat (Throat 2), migrate across the compact manifold M, and 
eventually annihilate the D3-branes. Here we work in the SUGRA limit to approximate the 
total migration time, in a two-throat geometry. As discussed in section the KPV/CDL 
decay times are so large that the migration times have little effect. This appendix, therefore 
serves largely to show that one may, in fact, ignore the migration time and to illuminate 
how the migration itself proceeds. 
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It will be assumed that the back-reaction of the migrating D3-branes is negligible, the 
proper velocity of the branes remains small, and that the majority of the travel time comes 
from the two throats (i.e. the time through the bulk of the CY may be ignored). 

We turn our attention once more to the metric © (with the overall scale of the manifold 
scaled back in). In 3(j] it was shown that the F 3 flux wrapped on the A-cycle smoothly 
deforms the tip of the conifold. Though we will find that the majority of the travel time 
comes from the tip of throat 1, we may ignore the most of the details coming from the 
deformation of the conifold since the motion is assumed to be radial. We will use the 
undeformed metric and, when working near the tip, multiply the warp factor h by an overall 
constant ~ 0.4 to account for the deformation 12 . This has little effect on the final result, 
however it was such a trivial correction it seemed silly not to include it. The warp factor, 
away from the tip, is 

h = ^Hr/r s ) (Al) 
. 27r(N + p) 1 

r - = r " ^^jp 1 - i> (A2) 

p = # of D3-branes ; r 2 = 3/2 5/3 . (A3) 

Due to the deformation of the conifold discussed above we will only be interested in the 
region, f = r s exp(l/4) < r < tq. Note that this avoids the naked singularity at r = r s . 
The action for the D3-branes is 



S 3 = -— / d 4 f y^d +/i 3 / C A ■ (A4) 

9s J JS(D3) 

The Ramond-Ramond potential C4 depends only on the radial distance r: 

C A = ^-dt A dx l A dx 2 A dx 3 (A5) 
9s 

where t = x°. Working in the gauge £° = r(t) and £ l = x l , the Lagrangian becomes, 

^ 3 l C = - — (y/t 2 -hr 2 ) + f(r)i. (A6) 
9s 

Assuming that the proper velocity is small, 

, lrH- 1 (h- 1 - g s f(r)) ■ 

u~ l C ~ i ysJ K JJ t A7 

2 g s 9s 



12 One finds this correction by comparing the "near tip" warp factor found in to the naive limit of the 
undeformed Klebanov-Tseytlin geometry. 
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It is easy to check that this a valid approximation. In particular, one only needs to consider 
throat 2, since this is where the D3-branes are moving fastest. One can show that p « 
(g s M 2 )/8 will insure that (|A"7|l is valid. 

Since t(r) is a cyclic variable, we know that dC/dt = —E/g s is constant, leaving us with 

i Pif 

E =2 { m )2 + V{r)] n^^ih-'-m) , (A8) 
and the travel time through a single throat is therefore, 

r dr 

At = ± - . (A9) 

J ro y/2(E-V(r)) 

The +/— corresponds to branes traveling into/out of the throat. 

Note that the time here is a coordinate time, but we will see that it is so small compared 
to decay times that we do not need to worry about conversion to proper time in the 4D 
Einstein frame. 



1. Throat 1 



The D3-branes, produced at rest in throat 1, feel a slight gravitational attraction to the 
D3-branes at the bottom of throat 2. Since the gravitational attraction is weakest while in 
throat 1, one suspects the majority of the migration time it comes from throat 1. We show 
below that, in fact, throat 2 can be ignored completely. This also justifies not including the 
travel time through the bulk of the CY. 

The Ramond-Ramond potential, C4 = h^ 1 / g s , is not effected by charge contained in 
throat 2. Physically, this is due to the charge being screened; mathematically, we are 
working on a compact manifold and may consider just the charge enclosed in throat 1. 
However, the geometry does, albeit slightly, know about what is happening in throat 2. 

In order to get an estimate of the migration time we will make the (perhaps bold) as- 
sumption that, as with most multi-pole solutions in gravity, the warp factor is changed by 
an additive factor, 

Mr)— /ii(r) + <y(r); S(r) = - ™™ ' (N 2 + p + 3g s M 2 /87c) . (A10) 

Here, the subscripts indicate the throat in which the fluxes are contained, p is the number 
of D3-branes in throat 2 (there are no D3-branes in throat 1). Expanding to first order in 
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Throat 1 




FIG. 3: Potential in throat 1 (Model 1) 

5, we see that 

277rr/ 

V(r) = - (2rQ _ (N 2 + p + 3fl s M 2 /87r) ^ 2 (r) . (All) 

Figure El shows the this potential. Note that the majority of the time will be spent at the 
tip of this throat. We have been unable to evaluate the integral (jA9|) explicitly. Numerical 
methods also proved difficult, due to the singular behavior of the integrand as r — > f. This 
results from the fact that the D3-branes are produced at rest. However, one may gain 
control of the situation by linearizing the potential near r = r and integrating away from 
the problematic singular point. It is in this limit that we multiply hi by the numerical factor 
~ 0.4 discussed above. Once a safe distance away from r = f, which for technical reasons 
coming from the linearization of (|A11|) we take to be r ~ f + zr , one may evaluate the rest 
of the integral ()A9j) numerically. For model 1, one finds that 

At! = 1.5 x 10 16 . (A12) 

2. Throat 2 

For completeness, we will determine time spent in throat 2, subsequently showing that 
it is of no importance. In order to deduce C4, recall that f -kdC^ ~ (N e ff —p). Plugging in 
the appropriate constants, this leaves us with 

r 5 h 2 d r f = (27W 2 )(A e// -p) . (A13) 
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Throat 2 




FIG. 4: Potential in throat 2 (Model 1) 

Using what we know from the KS geometry, we define / = h~ x + V, where the potential, 
V(r), satisfies 

_ (27*a»g.)p r*/L* 
T V (ln(r/r s+ )) " (AM) 

This can be integrated and gives a solution depending on exponential- integral functions. 
The potential is shown in figure EJ Numerically integrating (jA9|) for throat 2 gives 

At 2 « 9.9 , (A15) 

which is clearly negligible compared to (jA12j) . 

APPENDIX B: INCLUSION OF GRAVITY IN BUBBLE NUCLEATION 

In recent literature, there has been some confusion concerning the relation of two for- 
malisms for studying thin-wall instantons. The method of Coleman and De Luccia (CDL) 

nn 

|32J, |36( describes smooth instantons in the limit of large radius of curvature; this formalism 
was used by ((J to argue that bubbles will always nucleate in a dS background before the 
recurrence time. Alternately, however, we could imagine that the bubble wall is truly an 
infinitesimally thin membrane, such as a D-brane, with a delta function stress tensor. This 
type of configuration was studied by Brown and Teitelboim (BT) 12]. In a description 
of dS solutions in noncritical string theory jsj], uses the BT formalism to argue that 
there is a critical tension above which the bubble occupies more than half the original de 
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Sitter sphere and above which the decay time changes behavior as a function of the bubble 
tension. Additionally, claims that the decay time is of order the recurrence time at the 
critical tension. 

In this appendix, we show that the CDL and BT formalisms actually agree; this is 
reassuring, since even D-branes should be described as smooth objects in a complete version 
of string theory. Our results show that the decay time is always less than the recurrence time, 
as in [J 0, B, S| • Additionally, we confirm that more than half the original de Sitter sphere 
decays above the critical tension, but we show that (due to some technical considerations) 
the decay time is actually a smooth function of the bubble tension. We work in a 4D effective 
theory throughout. 

We begin by describing the two formalisms. In both CDL and BT, the nucleation/decay 
time is given by exponentiating the difference of the (Euclidean) bubble and background 
actions. Therefore, the exterior of the bubble, which is approximated by the background, 
contributes nothing in both formalisms. 

At that point, CDL note that, since both bubble and background have the same behavior 
at infinity, they can integrate some terms in the Ricci tensor by parts. After determining 
the bubble tension r as a functional of the potential, they find that the action is 



AS E = 2n 2 r 3 T + 
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where A± = n\V± are the potential outside and inside the bubble respectively (this is a 
combination of eqns (3.11) and (3.13) from 3^|). Minimizing this action with respect to the 
bubble curvature radius r gives the decay rate. 

On the other hand, BT cannot use the same integration by parts because the infinite stress 
of the bubble wall separates the interior and exterior regions. Instead, the bubble action 
must include extrinsic curvature terms; it is these terms that will explain the apparent 
contradiction between CDL and BT formalisms. The extrinsic curvatures of the interior 
and exterior regions are 



K 4 



— 3a± 



1 



3 



1/2 



vB2) 



where a± = 1 if the radius of curvature of the outside/inside region of the bubble is increasing 
toward the exterior of the bubble and is —1 if the radius is decreasing. Since the Ricci scalar 
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in the bubble is given by the cosmological constant, the action just becomes 

2tt 2 



AS r 277 2 r; ;! 4- ~r 6 (K- - K+) - \ (A_V_ - A+V 4 



(B3) 



The interior volumes for the bubble and background are given by (for either sign of the 
cosmological constant) 



V 4 




o"± ( 1 - 



3/2 



- 1 



1 



(B4) 



It is algebraically simple to see that the extrinsic curvature terms combine with the square 
root terms from the volume to give exactly the CDL action (|B1|) up to the signs a±. The 
reason found a different action is that they omitted the extrinsic curvature terms. 

Now we should see why the CDL result should actually have the signs a±. For a de Sitter 
background, the terms in square brackets of eqn (jBl|) come from integrals 



A 



±„2 



d£ r ( 1 —r 



(B5) 



which CDL evaluate by replacing d£ = dr(l — r 2 A±/3) 1 ^ 2 . However, as they note, 



y/3/A± sin[-y/A±/3£], so ^(r) is double- valued. In fact, the correct integral is 

dy y 2 



3 

A± 



(B6) 



'cr±(l-r 2 A±/3) 1 /2 

which just introduces a factor of a± in the (■ • • ) 3 / 2 terms. This precisely agrees with the BT 
results. In this paper, we will be concerned only with the case <r_ = 1, and a + = — 1 only 
for A + > and tension above critical. 

To find the radius of the bubble given the two cosmological constants and the tension, we 
could minimize the action with respect to r. However, it is easier to use the Israel matching 
condition across the bubble wall, which has trace K + — A_ = (3/2)k|t. The answer is given 
by and can be written as 



2 \ 2 
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AA = A_ - A_ 



(B7) 



It is tedious but straightforward to check that this matches the result from minimizing the 
action. The maximum radius occurs at critical tension 

1/2 

AM ) 
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and is 1/r 2 = A + /3 = l/R^s for positive initial cosmological constant. Note that for 
vanishing initial vacuum energy, gravity stabilizes the false vacuum for tension bigger than 
critical, as in (3^]. Also, for negative initial cosmological constant, the radius becomes 
infinite for tensions lower than critical. We will concern ourselves only with initial de Sitter 
spacetimes, so we do not face some of the concerns raised by [36] about decays of Minkowski 
and AdS spacetimes. 

We will finally write down the action for the bubbles: 



2ttV 




(B9) 



(BIO) 



While this is a mess, the reader should note that the sign of the last term is independent 
of the tension. That is because, for t < t c , a + = 1 but the quantity in the parentheses in 
the last term of (jB9|) is the square of a negative number, so the square root introduces a 
sign. For supercritical tension, that quantity is the square of a positive number, but then 
a + = —1. We have chosen to write the variables in this form in order to illuminate the 
dependence on the level spacing. Please see figure Elfor the qualitative features of the action 
AS as a function of A_, AA. While in some ways the physics depends more directly on the 
initial cosmological constant A + , in this paper we typically works with a fixed final A_, and 
AA depends on the same moduli that control the bubble tension. 

We should note that this action reduces to the known formulae in special cases. In 
particular, the result of CDL as quoted in KKLT, 

Sq 



AS? 



(1 + t c 2 /t 2 ) 



(Bll) 
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(a) Lines, from top to bottom, are at 
dL = -0.1, -0.3, -0.55, -0.75, -1,-1.5. 



(b) Dependence on both variables. 



FIG. 5: The bubble minus background action as a function of the cosmological constants. The 
variables are S = k\t 2 /S.S, L = A^/nfr 2 , and dL = AA/k|t 2 . 

is valid for all tensions when the final vacuum energy vanishes. A related result is that, for 
any A± > 0, as the bubble tension goes to infinity, the decay time goes to the recurrence 
time of the original dS. 



32, 361, that the final states are not 



As final comments, let us reemphasize, following 
maximally symmetric spacetimes but rather cosmological ones. In particular, decays with a 
negative final cosmological constant lead not to AdS but to a Big Crunch singularity within 
the bubble. Additionally, as mentioned in jj^j], these instantons are technically different 
from instanton decays of inflationary spacetimes. It seems reasonable that for sufficiently 
small decay rates that treating the initial spacetime as dS is a good approximation, but 
it remains an interesting problem to study decays of possibly more cosmologically relevant 
spacetimes. 
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